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The Motion of Viscous Liquid due to Uniform and Periodic 
Motion Maintained over a Segment of an Infinite Plane 
Boundary, 

By W. J. Harrison, M.A., Fellow of Clare College, Cambridge ; Assistant 

Lecturer in Mathematics, University of Liverpool. 

(Communicated by Sir J. Larmor, Sec. E.S. Eeceived October 14, — 

Eead November 21, 1912.) 

The author was led to a consideration of the problems whose solutions are 
presented in this paper by the attempt to solve a certain problem in the 
theory of lubrication. But it was thought that these solutions have sufficient 
intrinsic interest to warrant their separate publication as examples in 
theoretical hydrodynamics, the more physical applications to lubrication 
being reserved for further development. 

In the first three problems infinite viscous liquid is maintained in motion 
by uniform or periodic tangential motion of a part (a strip of infinite length) 
of its plane boundary. In the remaining two problems the liquid is further 
limited by a second plane boundary parallel to the first. 

A departure, which may be moi'e or less important, from the usual actual 
physical conditions is necessitated by the result that it is impossible to 
maintain the condition of no relative velocity at the boundary. The treat- 
ment of the problem with the condition of no slip leads to a discontinuity in 
the velocity of the liquid at the boundaries of the moving strip. But, on 
the other hand, from the nature of the problem this condition of no slip 
is an impossible physical state of affairs in the neighbourhood of two lines 
of the boundary. This complexity is avoided if we take the more general 
condition of slip with tangential traction proportional to the relative velocity, 
and then the case of no slip will be approached by taking the ratio of 
traction to slip to increase indefinitely. There are two other alternative 
suppositions which might be made and which seem to carry with them a 
more physical appearance. These will have to be examined before further 
applications are made of these solutions, but the solutions in their present 
form are given for their own interest. 

Problem 1. 

In this and the following calculations the motion is two-dimensional, and 
will therefore be considered to take place in the x, y plane ; accordingly, the 
third co-ordinate will never be mentioned. 
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Viscous liquid occupies the part of the x, y plane for which y is positive, 

being bounded by the axis of x. This boundary is fixed, except that between 

^ = ±$ it is given a uniform velocity th in its own direction. The effect 

may be produced by supposing the part of the axis — ^ <];:»<; ^ cut away, 

and underneath in contact with it a second boundary is placed which moves 

with the given velocity. It is required to determine the motion which is 

thus maintained. 

The motion is steady and is, moreover, assumed to be slow.* It follows 

/ 32 32 \ 2 
that the stream function satisfies the equation (5-2 + 5-2) '»/^ = 0. 

If we assume a disturbance depending on cos hx, then 

^^r = {Ae~^y + B?/e~^y) cos kx, 

where k is positive. This form of the stream function is appropriate to a 
velocity along the bounding plane varying as cos kx. Let the velocity of the 
boundary in its own direction be tjsu cos kx. 

It is impossible to satisfy tlie condition of no relative motion at the 
boundary owing to the discontinuity that would arise at the points 
^ = ± ?> 3/ = ill til® ^TLdl synthesis of the harmonic elements. Accordingly, 
the assumption is made that the boundary exerts a traction on the liquid 
proportional to the relative velocity. Thus 

-^—\ -7^ = Wk cos kx (y = 0), 

dy vy^ 

is the appropriate condition. If X be made very small the condition of no 
relative velocity is approximately satisfied in general. 

The remaining boundary condition is -^/r = (3/ r= 0). 

These lead to 

A=0, B = t«T;t/(l + 2X^). 

If the distribution of velocity of the boundary in its own direction be given 
by fix), then 

-I poo pCO /•/ \ 

Suppose /(a) = w, — ^ < ^ < ^, 

■^ The fact that the motion has to be assumed slow in order to obtain linear equations 
is a very serious hmitation on the interest of the solution ; but it must be noticed that 
when it is a question of motion between parallel planes which are very near together the 
non-linear terms are neglected by Osborne Reynolds in his paper on lubrication, for 
physical reasons, and the limitation of slow motion is thus removed. This consideration 
applies to the last two problems. 
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Hence 



^ = 2/™ 






y 



m 



tan 



-l^ + ^_Lfar.~l ?"~^ 



TT 



e¥ 



4- tan" 



1/ 



.00 







r — ^Yr [sin k (f + a?) + sin /<^ (^— a^)] f& 



A convenient method of evaluating the definite integral is to plot the 
integrand for particular numerical values of oo, y, ^, X. But it is to be 
noticed that if X be small, say X = 0*01, then, except near the boundary, 
i|r is given by the integrated terms to a high degree of approximation. 

In fig. 1 the stream lines are shown for the case X = O'Ol, f = 1. They 
are drawn for the equidistant values {\, f, 1, f, §)m/7r of the stream, 
function. They are, however, independent of X. 





X-8^ 



Fig. 2. 



i-zl 



Fig. 1. 



In fig. 2 the stream lines are drawn in the neighbourhood of ^ = 1, 3/ = 
for the same case. The method of evaluating the integral is similar to one 
to be mentioned later. The stream lines are drawn for 

yjr = (0-03, 0*06, 0-12, 0-24) mjir. 



How 



Distribution of Velocity along the Axis of x. 






wf sin>fc(^4-^) + sin/^(^— a?) 



^(1 + 2X7^) 
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/•QO 

Jo 



sin ah 



dk 



.00 

ah 



sin(m— a5) 



m 



dm (provided a > 0) 



Si {ah) ) cos {ah) + Ci {ab) sin {ah). 
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Hence 



m 
m — — 

TT 



/|"_Si(9)cos(9 + Ci(9sm^ 



TT 

,2' 



+ ( ^ — Si ^ I cos cj> + Gi (f) sin q^ 



0<x<^, 



m 

TT 



A. 



I ^ Si ^j cos ^ + Ci sin ^-f ^-Si <!>) cos ^-~Ci ^ sin ^ 

where ^ = (^'j-x)/2\ (j) = | ^—x\/2\. 

The distribution of the velocity along the boundary, of which the moving 
section is bounded by a? =+ 1, is exhibited in the following table. The 
evaluation has been performed by the aid of Glaisher's Tables of the Sine 
and Cosine Integrals. 



I 

A. 


0? = 0. 


0-6. 


1. 


1-4. 


2. 


3. 


0-01 
0-1 

1 


0-987 

0*88 

0-452 


0-98 

0-833 

0-416 


0-497 
0-469 
0-303 


0-013 
0-101 
0-183 


0-039 
0-119 


0-017 



The table shows the ratio of the velocity to m. 

In addition to the above values it is interesting to note, in the case of 
X = 0*01, how rapidly the velocity changes in the neighbourhood of ^ = 1. 

^= 0-9 u = 0-937 GT, 

X = 0-99 II = 0-713 w, 

x = l -2^ = 0*497 t«T. 

From the distribution of the relative velocity of the sliding piece and 
the liquid in contact with it, the traction {/jud^^/df, 2/ = 0) exerted by the 
one on the other can be directly obtained, 












dy'r 



y = o, 



where fju is the coefficient of viscosity. 

The total traction exerted by the liquid on the sliding piece per unit of 
its length is given as follows : — For the values 0*01, 0*1, 1 of \, traction//>6T^ 
has the values 31, 1*6, 0*56, the breadth being two units. 

The total rate of flow in the direction of motion of the sliding piece is 
2m^/7r, and this is independent of X. This is to be explained by the fact 
that, although the sliding piece exerts less traction on the liquid the greater 
the value of X, the liquid itself experiences less resistance from the remainder 
of the boundary. 
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Problem 2. 

From the foregoing solution it is possible to derive immediately the 
stream lines where there are two such sliding pieces in different parts of the 
boundary. 

In fig. 3 the stream lines are drawn for the case when the part of the 
boundary 1 <Cx <^Z is given a velocity m in the direction of x, and the 
part — l>a:5>— Sis given an equal velocity in the opposite direction. The 
stream lines are drawn for the equidistant values (0*5, 0*75, 1, 1*25, 1*5) x^jir 
of the stream function. 





Fig. 3. 



Fig. 4. 



In fig. 4 the stream lines are drawn for the case in which the velocitf 
of the second sliding piece is half that of the first. The lines are shown 
for the same values of '^ as in the previous case, except that ^ = 0*25 mjw 
is drawn as well. 

ProUem 3. 

We proceed to consider the case in which the sliding piece, instead of 
having a uniform velocity, is given a periodic velocity mQo^pt, 
The stream function must now satisfy the equation 

where v is the kinematic coefficient of viscosity. 
The appropriate solution is 

'\^ = (Aa~'^^ + Be"^"*^) cos hxe^P*, where m^ = k^+ip/v, 

leading to m = q-\-i/j., 

VOL. LXXXVIIL — A. 
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q being taken positive. Also, since we are concerned with a maintained 
motion, p is real. 

The axis of x is a stream line, and therefore 

A + B = 0. 

This solution is appropriate to the harmonic distribution of velocity 
GT/fc cos kx e^^^ of the boundary. Hence the remaining boundary condition is 

^-^X^ == wk cos kxe^P^ (3/ = 0). 

This leads to — (^ + XP) A— (m + Xm^)B = sr/fc. 

Hence A = —B = mk/iin—k + iXp/v), 

Taking real parts, we find 

yfr = tij^F (k) cos kx, 
where 

F(^) = [e~^^ (P^in pt-^Q cos pt)--'e~^^ {V {G0HfjL2/^mpt — ^m fjuy cos pt) 

+ Q (cos fjby cospt -f sin fiy sin^^)}]/(P^ -f Q^), 

where P = fi + Xp/v, Q = q—k. 

The solution, following the notation of Problem 1, is 

'Jr = - I [m\k(^+x) + smk(^---x)]---Y-^ dk 

The work of calculating the stream lines for particular values of the 
constants involved is a matter of some complexity. I have performed it for 
one case, but not sufficiently to give more than a rough idea of their form. 
But it was possible to see that a vortex is periodically formed and destroyed 
in the neighbourhood of the sliding piece in this particular case. It would 
be interesting to trace the stream lines for a range of values of the time, and 
so exhibit the growth and decay of this vortex. But at the present time 
I have not the requisite leisure. 

Distribution of Velocity over the Sliding Piece.— it can easily be shown that 



f-g^)^^^ = Jj [sin^(? + ^) + sin^ T F+QV^^^^^ 



+ ^~^2«i^i^^ 



dk 
~k 



The numerical evaluation of this integral presents feiatures of some 
interest. Apart from x the variables are ^, X, plv * ^ has been taken equal 
to one unit as before, X = 0*01, and the integral has been evaluated for 

^Iv = 10~i, 1, 10, 102, 10^ 10^ 10«. 
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The method of evaluation depends upon the magnitude of p/p, and one or 
two illustrations of the methods employed will be given. 

(a), p/v == 0*1. 



r P . ,,^ ^dk f2 P . ,,^ ^dk rsmk(P-{-x)dk 



k ' 



since Q is very small for ^>2. At the same time /a =:j?/2^i/, so that the 
second integral is capable of evaluation by the aid of the tables of sine and 
cosine integrals as in similar integrals in Problem 1. The first integral can be 
evaluated by plotting the integrand. 

In this same case the coefScient of sin pt can be evaluated by plotting the 
integrand, since it decreases rapidly with h 

(h). p/p == 10. 

The coefficient of m cospt (for ^ = 0) involves an integral which can be 
replaced by 

'(0-223 -f 0105 Z: + 0-0063 F) ^ dk + f^^^^^ 
10^ Jg bO-i-k 

to a good degree of accuracy. These integrals are easily evaluated. 
(c), p/p = 10*. 

The coefficient of m co^pf (for x = 0) involves an integral which can be 
replaced by 



j' 

Jo 






100 „,•„ h 

(0-005 + 0-0000184 k) ^ dk 





k 



Or (d). p/v ^ 10*. 



^.)eosi.^+^,sin^^_ 



for any value of x between ±^, except near these limits. 
When p/p is very great ( ^ ) =0. 

Eepresenting the velocity at the centre of the sliding piece by 
X cosp^4- Y sinp^, we have the following values for X and Y. 



I>l^ 


10-1. 


1. 


10. 


101 


103. 


101 


101 


X ' 

Y 


0-99 
0-001 


-99 

0-005 


-97 
0-025 


-93 

0-06 


0*79 
0-16 


0-5 
0-21 


0-07 
0-06 





X decreases from 1 to 0, Y increases from to a maximum value and then 
decreases again to zero. The velocity at the point ic = I f is 1 or 2 per cent. 

2 
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below that at the centre, and at x = ^,X. and Y have practically half their 
values at the centre. 

The maximum (as regards the time) integral traction on the sliding piece 
per unit of its length can easily be found. It has the following values in 
units of the system of measurements employed. 



J^l^ 


10-1. 


1. 


10. 
14 


10-^ 


W\ 


10^ 


10^. 

i 


Traction 

fJLVT 


3 


3 


1 

23 1 

1 


55 


114 


190 



For water, p/v = 10^ corresponds to a period of 0*6 sec. in C.G.S. units, 
so that for this period the traction is about 18 times that for steady motion. 

Problem 4. 

This is the same as the first problem, except that the liquid is further 
bounded hj y = b, which is fixed. The corresponding boundary condition is 
taken to be that no relative motion is possible over it. 

A possible stream function is given by 

A^r = { (A + Gy) sinh % + (B + T)y) cosh ky} cos hx. 
The conditions to be satisfied are 



^jr = 0, ■^"~'^52 ~ ^^ ^^^ ^^^ (^ 



0); 



t = o, 

These lead to B = 0, and 
A D 



1 = (. = .). 



Icb^ 



kh — I sinh 2 kb ^inh^kb 



VSfc 



mic 



- k%^ - 2 mx + k\ sinh 2kb-i- sinh^ kb ^ (k) 
Using; Fourier's theorem, we obtain the stream function in the form 



.|r= ~J [sinZ^(^+^) + siny^(^~-^)]^-^ , 



where 



^ (^) = [ — kb^ + (A:& — J sinh 2kb) y] sinh ky + y sinh^ kb cosh ky. 

In fig. 5 the stream lines are drawn for the case of ^ = 1, & == 2, X =1. 

Since calculating ^fr for the unpractical case of X = 1, I have lighted upon 
a method of evaluating the integral for X = O'Ol, but the advantage of 
recalculation is not great in the case of b so great as 2. 
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Fig. 5. 



Velocity of the Liquid at the Centre of the Sliding Piece, — The expression for 
the velocity of the liquid along the boundary 3/ = is 






.00 



— [sin h (^ +^) 4- sin h (^— ^)] dh 







h(f>{h) 



The ratio of the velocity at the centre of the sliding piece to cy is given 
below : — 



A. 


5 = 1. 


0-1. 


0-01. 


0-01 
0-1 


0-96 
0-71 


0-72 
0-2 


0-2 
0-025 



When h is small, it is easily shown that the velocity at any point of the 
sliding piece is hwl(b-\-4X), except that the approximation breaks down when 
^—x is small, and that sit x = ^ the velocity drops to half that at the centre. 

Hence the traction exerted by the liquid on the sliding piece is (per unit 
length) 



2 



tjy< 



\ 



& + 4X/ 



M 



b + 4:\' 

As h is decreased indefinitely, this approaches the limit 2//,ct^/X. This is 
the same limit as in the case of the rapid periodic motion of the last case. 

Frohlem 5. 

The liquid is bounded by two parallel planes as in the previous case, but 
the sliding piece is given a periodic velocity m aospt. 
The appropriate typical harmonic solution is 

i|r = ( A cosh ky + B sinh ky + G cosh my + D sinh my) cosh kx e*^*, 

with the same notation as in Problem 3. 
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The solution follows as in previous cases, and I proceed at once to discuss 
the velocity of the liquid at points on the sliding piece. 

Consider '~j~'fQ<^) ^T^- 

V j 

If /(/c) vary slowly with h within the range < ^ <[ 20 (say), the value of 
the integral is -|7r/(0) to a very fair degree of accuracy. 

It is assumed that /(A;) is decreasing and one-signed. The rate of variation 
I need not be very slow for the approximation to hold, for example / (k) may 
decrease practically to zero in the interval. The approximations which 
follow are based on this consideration, but it is not claimed that the inter- 
mediate results in the range are more than roughly correct. 

where f{h) is the real part of e^^* divided by 

1 _L ^^^ ''^^ ^^^^-^^ ^^ cosh onh—h cosh kh sinh mh 

V ' 2 km{l — cosh kh cosh mb) + (P + m^) sinh kh sinh mh ' 

Now, provided h be small, the variation in f(k) is small for moderate 
values of k. Accordingly /(/^) is replaced by/(0). 

Now, when k = 0, m = q + i/ji, where q — fx — (jpj'lvf — ol (say). 

Put a6 " 7], \0L = ^. 

It is found that 



-^ 



where 



w (E cos jpif + S smpt)/T, 



E = 2 — 4 cosh 7] cos ?; — 2 ?;f sinh ?; sin ?; + cosh 2 ?; + cos 2 t] 

— (2?;--^) (cosh 7] sin?; --sinh?; cos ?;) + J(2?; — ^4"^^^ sin 2?; 

— J (27] — ^-—7]^^^) sinh 2?; + !?; (?? — ^) (cosh 2?; — cos 2?;), 

S = ^[?;(cosh 2?; + cos 2?; — 2 cosh?; cos ?;) — |(l-f ?;^) sinh 2?; 

— I (1 —?;^) sin 2?; + (cosh ?; sin ?; + sinh 7] cos ?;)], 

T = 2— 4 cosh?; cos ?;— 4?;^ sinh?; sin?;— ■2(?; — ^) (cosh ?; sin ?; 

— sinh?; cos ?;)+ (?; — ^) (l-f-^?) ^i^ 2?; — (?;~-^) (l — ?;^) sinh 2?; 
4- (1 + vV) (cosh 2 ?; + cos 2 ?;) + i (?; - ^f (cosh 2 ?; - cos 2 ?;). 

(a) Approximation when och is small.-— 

The velocity at the centre of the sliding piece is equal to - — - — vs cos pt, 

-{-4: Ay 

the coefficient being the same as in the case of steady motion. 

This approximation is valid if (1) h = 0*1 cm., a ^ 1, which gives for water 
a period of 5 mins. and over, (2) h = 0*01 cm., a < 10, which gives for water 
a period of 3 sec. and over. 
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(b) Approximation when ocb ^ 4, so that cosh ccb and sinh ab can be neglected, 
compared with cosh^ ccb, etc. — ■. 

This applies in the case of water to periods of 0*2 second and under, and 
0*002 second and under, for the cases & = 0*1 and b = 0*01 respectively. 
The velocity at the centre of the sliding piece is 

The following table includes the w^hole range of values of ocb, for \ = 0*01 
b = 0*01, in which case tj = ^, and the above formula is simplified. 



ab 


0-01. 


0-1. 


1. 


2. 


3. 


4. 


5. 


6. 


10. 


R/T 

-S/T 


0-2 



0*2 



0-199 
0-011 


0-19 

0-088 


0-157 
0-082 


0-121 
0-07 


0-097 
0-064 


0-083 
056 


0*05 
0-05 



It follows that the traction exerted on the sliding piece bears a ratio to 
that exerted in the steady motion which lies between 1 and (b + 4: X)/4 \. 

It will be noticed in Problem 3 when the period of the motion is very 
small, and in both of the last two problems when b is small, that the 
Condition of no slip at the boundary is not approximately satisfied ' for 
\ = O'Ol. If \ be taken sufficiently small, then it is easily shown that the 
maximum total traction per unit length on the moving piece has a value 
2/iim^{p/v)i, in the case of very rapid motion in Problem 3. In the steady 
motion of Problem 4, when the distance b between the two planes is 
very small, the integral traction has the value 8/jum^/b, In the periodic 
motion of Problem 5, if b be small and b(p/2v)^ < O'l, the integral traction 
depends upon b and not upon p, and has the maxiinum value S/jlvs^^/I). If, 
however, b be small and p sufficiently great, the traction is independent of h 
and has the value 2fim^(p/v)i. This can be shown from the approximate 
expressions for E/T, S/T, which are given above. 

Thus in the case of periodic motion between two near parallel planes the 
integral traction is equal to or greater than Sfjum^/b, and less than or equal 
to 2fjim^{ylv)^, provided this second limit is greater than the first. 



